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Possible Ac Ac hadronic molecule is investigated in the one-pion-exchange potential model. In the 
study with this model, the heavier meson exchange elfects are encoded into a phenomenological 
cutoff parameter and couplings to the nearby EcEc, ScS^, and E'Ec channels are essential. From 
the numerical results, we find that a molecular bound state of two Ac's is possible, where the tensor 
force plays a crucial role, although the binding energies are sensitive to the cutoff parameter. 

PACS numbers: 12.39.Pn, 12.40.Yx, 13.75.-n, 14.20.Lq 

Hadronic molecules are loosely bound states of hadrons, whose inter-hadron distances are larger than the quark 
. confinement size. The deuteron is the well-established molecule composed of a proton and a neutron. The triton, 
' hypertriton and so on are also regarded as molecular bound states of the light baryons. Recently observed near- 
(— I , threshold charmonium-likc mesons triggered lots of studies on the molecule problem in the heavy quark realm. In 
' this note, we consider the molecule problem of two Ac's. 

Compared with light baryons, the heavy quark baryons are more likely to be bound. One reason is the larger 
reduced mass of the system. The relatively small kinetic term in the Hamiltonian is advantageous for the bound 
state. The other reason is the heavy quark spin symmetry and thus the importance of channel coupling. In the limit 
■ of infinitely heavy quark, QCD interaction manifests heavy quark flavor symmetry and heavy quark spin symmetry. 
' The latter symmetry leads to degenerate Sc and S*. In the real world, the mass difference between Ec and E* is 
O |. indeed smaller than that between E and E*. Since the coupled channel effects become important when two channels 
JL^ ' are closer, it is necessary to include such effects in the study of heavy quark baryon interactions. We explore the 
^ , importance of such effects in the AcA^. molecule problem. 
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Channels 1 2 3 4 5 

J"^=0+ AcAc('So) ScSc('go) KK{'So) KKCDq) ^cKCDqJ 

TABLE I: The S'-wave AcAc state and the channels which couple to it. 
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I The quantum numbers of the 5-wave AcAc arc I{J^) = 0(0+). Here we consider in total five channels which are 
■ given in Table HI The wave function of the 5-th channel is taken to be 

^ ■ 1 / 

o : iece:) = {[^cK]s==, - p:ec]^5=o) , (i) 



where the minus sign comes from the exchange of two Fermions. 

In Ref. [l|, we have explored AcA^ system by using both one-pion-exchange potential (OPEP) model and one-boson- 
exchange potential model where exchanges of scalar and vector mesons arc also included. It is observed that the two 
^ . models may give consistent binding energies and corresponding radii. In the OPEP model, the contributions from 
shorter distance interactions are encoded in a phenomenological cutoff parameter. Since there is no pion exchange in 
the AcAc channel, the possible binding solution must result from the coupled channel effects. In the present study, 
we use OPEP model and investigate whether long range interaction may lead to a molecular bound AcAc state. 

Besides the five channels in TableHl '^ccN may also contribute. But its contribution may be important only at short 
distance since the exchanged mesons between 'E^ccN and any channel in Table |T] are much heavier. Here we neglect 
the 'E.ccN channel as we arc considering the possibility of loosely bound molecule. It was proposed that a bound state 
may exist also in the 'E.ccN system in Ref. (2|. 
The interaction Lagrangian reads Q 

C^nt = gitr{BQ-i^,-i^A''Bfi) + [.g2tr(B67M75^''S3) + h.c] + [53tr(B* A^Be) + h.c] 



[qMbTA^B-^) + h.c] + gMBVlt^l^A^^Bl.) + g^B-^^^^^A^^ B-^), (2) 
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Bq is similar to Bg, and = ^[^^(^^iO + with ^ = exp[^] is the axial vector current. The decay constant 



in the chiral limit has the value / = 92.3 MeV. In the heavy quark limit, the heavy quark spin symmetry requires 

a few relations for the coupling constants, i.e., 33 = -^51, 35 = ~§5ij 54 = ~VSg2, ge = 0. If one further uses the 

quark model symmetry, one has gi = 52- The relative phase between gi and g2 is actually irrelevant From 

the decay widths of Sc and S* we obtain \g2\ = 0.598 and \g4\ = 0.999 after averaging over the different charge 
states. In the numerical evaluation, we consistently use the following values and phases 



32 = -0.598, 54 0.999, gi 
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With the above Lagrangian, one derives the non-relativistic potentials. To incorporate the extended structure of 



baryons, a monopole type form factor F{q) 



'AT- 



is introduced phenomenologically at each interaction vertex 



where q is the pion 4-momentum. In principle, the cutoffs at the vertices AcSc^r, AcS*7r, Sc^cTt, Yic'S*tt, and I]*I]*7r 
are different. To reduce the number of parameters and simplify the calculation, we use the approximation that these 
five cutoffs are equal and we label this common cutoff A^. A^ is poorly known but its value, around 1 GeV, may be 
comparable to the nuclear models [H, Q . This parameter also plays a role of compensation for the short and inter- 
mediate range interactions. In our study, we treat it as a free parameter and discuss whether the AcAc molecule-like 
bound state is possible within the reasonable domain of A^. We may denote the potential in the following form, 



Oi ■ 02Yiir) + OtenHsir) 



(5) 



where i, j from 1 to 5 are the labels of the channels, C{i,j) is the coefficient containing the coupling constants, Oi 
(O2) is the Pauli matrix ct, the transition spin matrix St (or its Hermitian conjugation) explained below, or the matrix 

— (Oi • O2) is the tensor operator, and Yi, and relevant functions are defined 
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The transition spin for the Rarita-Schwinger field m'^ is defined through = S^^, where $ denotes the spin 
wave function of E* defined by 



$(3/2) = (1,0,0,0)^, $(I/2) = (0,1,0,0)^, $(-I/2) = (0,0,l,0)^, $(-3/2) = (0, 0, 0, 1)^ 



Explicitly, the time component of vanishes in the static limit and the other components are 
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In deriving the potentials, we have neglected the (5-functional term of the central force since we are considering 
molecule-like bound state problem. In the above potentials, m (A) is not always (A^r). In the transition potentials 
^15 J 1^25, Vijs, and V45, non- vanishing time component qq of the pion 4-momentum may be a better approximation. 
In these cases, we have m = ^/rn^'^~q^ and A — ^y A^ — q^. For the value of \qo\, we use (ms. — for V15 and 

V25, and (m|;. — ml )/(4ms*) for V35 and V45. Note that there are two terms in the final potential V55 due to the 
antisymmetrization given in Eq. ([l} , 
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We use go = in the gigs part and \qo\ = ms* — in the IgsP part. The above Idol's are derived in the static hmit 
of the heavier side, initial states or final states. 

For the hadron masses, we use = 137.27 MeV, tua, = 2286.46 MeV, ms, = 2453.56 MeV, and = 2517.97 
MeV In Fig. [T](a), (b), and (c), we plot diagonal and transition potentials of 5- wave case, S-D transition case, 
and D-wave case with the cutoff parameter = 1.0 GeV, respectively. From the diagrams, it is obvious that the 
tensor forces are strong and thus the coupled channel effects may be important. Another observation is that all 
the diagonal potentials are repulsive in this simple model. Therefore, the binding solution would result purely from 
coupled channel effects. 






275.0- 




220.0- 






> 


165.0- 


CD 






110.0- 








55.0 




0.0- 




-55.0 



FIG. 1: Diagrams (a), (b), and (c) show diagonal and transition potentials of S-wave case, S-D transition case, and _D-wave 
case with the cutoff parameter — 1.0 GeV, respectively, (ij) denotes the potential Vij(A.^, r). The last diagram shows wave 
functions with the cutoff parameter = 1.1 GeV. 



A^ (GeV) 


1.0 1.1 1.2 


1.3 


1.4 


1.5 


B.E. (MeV) 


3.39 14.45 35.44 


68.37 


115.06 


177.07 


VW) (fm) 


2.0 1.2 0.9 


0.7 


0.6 


0.5 


Prob. (%) 


(97.4/0.2/0.2/ (94.3/0.5/0.5/ (90.7/1.1/1.0/ ( 


86.8/1.8/1.8 


/ (82.8/2.6/2.8 


/ (79.0/3.4/3.9/ 




0.6/1.6) 1.3/3.4) 2.0/5.2) 


2.6/7.0) 


3.3/8.5) 


3.8/9.9) 


D-wave prob. 


2.2% 4.7% 7.2% 


9.6% 


11.8% 


13.7% 



TABLE II: Binding solutions for the coupled Ac Ac system with 5-channel contributions. Binding energy (B.E.) is given with 
relative to AcAc threshold. The probabilities correspond to AcAc(^S'o), EcEc(^So), T,*T,*CSo), T,*T,*{^Do), and EcE*(^_Do), 
respectively. 

It is not difficult to solve the coupled channel equations using the variational method ■ We obtain the numerical 
results in Table HIl with all the five channel contributions. If we drop the 13- wave channels, we do not find any binding 
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(GoV) 


1.0 1.1 1.2 1.3 1.4 1.5 


B.E. (MeV) 
(fm) 
Prob. (%) 


0.07 3.86 15.06 35.90 68.39 114.25 
11.3 1.9 1.1 0.8 0.7 0.6 
(99.6/0.0/0.0/ (97.0/0.5/0.1/ (93.4/1.5/0.2/ (89.0/3.1/0.4/ (84.2/5.2/0.6/ (79.3/7.6/0.9 
0.4) 2.4) 4.9) 7.5) 10.0) 12.2) 



TABLE III: Binding solutions for the coupled A^Ac system with only one D-wave channel. Binding energy (B.E.) is given with 
relative to AcAc threshold. The probabilities correspond to AcAc(^S'o), EcEc(^S'o), E*E*(^S'o), and EcE*(^_Do), respectively. 



solution. To see the importance of the tensor force, we include the contribution of only one D-wave channel. For 
the case of E*i;*(^£)o), one does not find binding solutions for < 1.3 MeV. The results for the case of Y,c^*{^Dq) 
(but without S*S*(^Do)) are given in Table Hill It is clear that the tensor force from S-D wave mixing is essential 
in getting binding solutions. From these results, one concludes that the channel I]cS*(^Do) pl^iys a more important 
role than Y.*j:i{^Do). 

As an example of the five channel solutions, we show the wave functions with A^r — 1.1 GeV in Fig. [T] (d). To see 
the sensitivity of the binding energy on the cutoff A^, we present a diagrammatic form for the results from Tabic HIl 
in Fig. [21 where we also show the binding energies for the uncoupled channel (with S-D mixing but without 

mixing or AcAj, channels). There is no binding solution in other uncoupled channels since the potentials are all 

repulsive. From Fig. [51 if A^r > 1.0 GeV is reasonable, a bound state is possible although the diagonal potentials are 
all repulsive and there are no binding solutions in individual channels. This indicates the importance of the tensor 
force. 

Although we do not have enough information to determine the cutoff parameter for the heavy quark baryons, we 
have interesting results in a reasonable range of Ajr. From the experience of nuclear force, the cutoff should be around 
1.0 GeV or larger, depending on the model. For the heavier hadrons, the extended structure is smaller and the 
cutoff parameter should be accordingly larger. In this study, the solutions corresponding to A^r = 1.0 ~ 1.2 GeV are 
molecule-like because the bound state is not so deep and the inter-hadron distance is not so small. A larger cutoff 
results in a tightly bound state and the OPEP model may be inapplicable any more. In all, it is possible to have a 
bound state of two Ac's while the binding energy is not determined precisely with the present approach. Wc hope 
that future studies may specify the binding energy of such a molecule state. On the experimental side, finding the 
double-charm AcAc bound state will be a challenging subject at GSI, J-PARC, RHIC, or Belle. 

In short summary, we have investigated the S-wave AcAc molecule problem by including the coupled channel effects 
caused by Ec and E* in a one-pion-exchange potential model. The couplings to the Z?-wave channels EcE* and S*!]* 
are crucial in binding two Ac's. The results are sensitive to the cutoff parameter A^. If the model cutoff around 
Att = 1.0 ~ 1.2 GeV is reasonable, one gets a molecule-like solution. 

Note: In a recent paper Q, the authors also studied AcAc system but they did not consider the excited S* 
contributions. The omission of D-wave channels results in different conclusions. 
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FIG. 2: Sensitivity of the binding energy on the cutoff parameter A^. 



